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1. Introduction 

Despite the great advances made over the last sixty years in the study of integrable 
quantum models, evaluation of their correlation functions still poses a formidable problem. 
Quite recently there has been significant progress in this direction: the group at RIMS 
succeeded in deriving integral representations for some correlation functions of the Heisenberg 
XXZ model defined by the hamiltonian (1.1) for A > 1 by taking advantage of the infinite 
quantum affine symmetry of the model on the infinite chain'^'^l . The isotropic (XXX) limit 
A — > 1 was obtained in [3,4]. These integral representations are most powerful for studying 
the short distance asymptotics of correlators, whereas it is not obvious how to extract the 
long-distance behaviour. Also it is not straightforward to extend this approach to the critical 
regime — 1 < A < 1 or to include an external magnetic field. 

Precisely these issues can be very naturally addressed in the framework of a different approach 
to studying correlation functions in integrable models, which was carried out in [5-13] for 
the example of the Bose gas^. We call this method the Dual Field Approach (DFA). The 
DFA is directly based on the (algebraic) Bethe-Ansatz solution of the model and thus is 
applicable to a large variety of correlation functions and integrable models. It allows to 
derive explicit expressions for the large distance asymptotics of correlation functions (even 
at finite temperature) , and the inclusion of an external magnetic field poses no problem. The 
DFA thus nicely complements the approach of the RIMS group. In a series of papers we will 
apply the DFA to the Heisenberg XXZ and XXX chains at zero temperature in a magnetic 
field h, i.e. the hamiltonian 



There are four main steps in the DFA: First the model needs to be ''solved" by means 
of the Algebraic Bethe Ansatz. Then one uses this solution to express correlation functions 
in terms of determinants of Fredholm integral operators. In step three these determinants 
are embedded in systems of integrable integro-difference equations (IDE). Finally the large- 
distance asymptotics of the correlators is extracted from a Riemann-Hilbert problem for 
the IDE's. As the computations for the various steps are rather involved we will only deal 
with the first two steps here, i.e. review the known Bethe Ansatz solution for the XXZ and 
XXX chains and then derive determinant representations for correlation functions. In two 
following publications we will present steps three'^^^ and four. 



A detailed and complete exhibition of this work can be found in the book [14] 
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2. A Short Review of Algebraic Bethe Ansatz 



Let us review a few main features of the Algebraic Bethe Ansatz (ABA) for both XXZ 
and XXX Heisenberg magnets, in order to fix notations for things to come. The XXX case 
can of course be obtained by taking a certain hmit of the XXZ case, but in practice this is 
more difficult than treating the XXX case separately from the beginning. Thus we will treat 
both cases on an equal footing throughout this paper. 

Starting point and central object of the Quantum Inverse Scattering Method is the R- 
matrix, which is a solution of the Yang-Baxter equation. For the case of the XXZ and XXX 
models it is of the form 



i?(A,/x) 



//(//, A) \ 

g{ii,X) 1 

1 ^(/i,A) 

\ fill, X) J 



(2.1) R 



where for XXZ 



sinh(A — u + 2ir7) , isin(2r?) 

/ A, /X = . " , ^ A, = . " . , 2.2 f gxxz 

smh(A — /i) smh(A — /xj 

and for the XXX-case 

/(A, /.) = ! + , y(A, li) = . (2.3) f gxxx 

A — jJ, A — jl 

The R-matrix is a linear operator on the tensor product of two two-dimensional linear spaces: 
R{lJi) e End{<C^ (g) (D^). Prom the R-matrix (2.1) one can construct an L-operator of a 
"fundamental spin model" I^^l by considering the matrix R{fj,)]l, where 11 is the permutation 
matrix on (D^ (D^, and then making it into an operator- valued matrix by identifying one 
of the linear spaces with the two-dimensional Hilbert space 7i„ of SU (2)-spins over the n'th 
site of a lattice of length L 

The Yang-Baxter equation for R implies the following intertwining relations for the L- 
operator 

R{X - /x) (LJA) ® LM) = (LM ® KW) - ■ (2-5) intL 

Prom the ultralocal L-operator the monodromy matrix is constructed as 

T{fi)^LML,_M...LM-(^ctJ^ ^1^)) ■ (2.6) monodr 

The intertwiner (2.5) can be lifted to the level of the monodromy matrix 

R(X - At) (r(A) ® r(//)) = (T(ii) (g) r(A)) R(X - n) . (2.7) intl 
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Below we will repeatedly use especially the following matrix elements of (2.7) 

[i?(A),5(/i)] = 0=[C(A),C(^)] 
[5(A), C{^^)] = giX, /x) (D(A)A(/i) - Z^(/x)A(A)) 
D{i2)B{X) = /(A, i2)B{X)D{i2) + g{i2, A)fi(^)D(A) 
A{i,)B{X) = fifi, X)B{X)A{i,) + g{X, fi)B{i,)A{X). 



(2.8) intTme 



By tracing (2.7) over the matrix space one then finds that the transfer matrices t{ii) — 
tr{T{ij)) — A{ijl)+D{ii) commute for any values of spectral parameter fx, i.e. [t{ijl),t{u)] — 0. 
From this it follows that the transfer matrix is the generating functional of an infinite number 
of mutually commuting conserved quantum operators (via expansion in powers of spectral 
parameter) . One of these operators is the hamiltonian 



d 

H^^, = -2zsm(2?7) — ln(r^^^(/i)) 



- 2Lcos(2r/) - 2hS' . (2.9) Hxxz 



Below we also make use of some properties of inhomogenous XXX and XXZ models, which 
are constructed in the following way: we first note that the intertwiner for the L-operator 
(2.5) still holds, if we shift both spectral parameters A and /i by an arbitrary amount u^, i.e. 

R{X - fi) (LJA - i/J ® L^ifi -iyj) = (LJ/. - i/J ® LJA - i/J) R{X - fi) . (2.10) intL2 

The reason for this fact is of course that the i?-matrix only depends on the difference of 
spectral parameters. We now can construct a monodromy matrix as 

T^nkW = L,{X - ^jL^_i(A - ...L,{X-v,)^ (^^^^j ) . (2.11) Tinhorn 

The ABA deals with the construction of simultaneous eigenstates of the transfer matrix 
and the hamiltonian. Starting point is the choice of a reference state, which is a trivial 
eigenstate of r(/i). In our case we make the choice |0) = | tit •••!) = ®n=il T)n> '^^ 
choose the completely ferromagnetic state. The action of the L-operator (2.4) on | t)„ can be 
easily computed and implies the following actions of the matrix elements of the monodromy 
matrix for the XXZ case 



A{fi)\0) = a(/i) |0) , a(/i) = (sinh(yU — ir]))^ 
D{ii)\0) = rf(/i) |0) , d{ix) = (sinh(/x + ir]))^ . 



C(/i)|0) = 0, 
B{fi)\0)^0, 

whereas in the XXX case 



(2.12) ADxxz 



a(/^)=(/^-^j , = (/^+ ^ j ■ (2.13) adxxx 
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Prom (2.12) it follows that -B(A) plays the role of a creation operator, i.e. a Fock space of 
states can be constructed as 



N 



*jv(Ai,...,AJ = n^(A,)|0) . (2.14) states 

i=i 

The requirement that the states (2.14) ought to be eigenstates of the transfer matrix r(//.) 
puts constraints on the allowed values of the parameters A„: the set {A^ } must be a solution 
of the following system of coupled algebraic equations, called Bethe equations 

^(Aj) fLi/(Aj,AJ 

These equations are the basis for studying ground state, excitation spectrum and thermody- 
namics of Bethe Ansatz solvable models. For the case of the XXZ model with A > — 1 (the 
case we are interested in here) it was proved by C.N. Yang and CP. Yang in [16,17] that 
the ground state is characterized by a set of real Xj subject to the Bethe equations (2.15). 
Without an external magnetic field {h = 0) their number is N = L/2. In the thermody- 
namic limit the ground state is described by means of an integral equation for the density 
of spectral parameters p(A) 



2np{X)- f dpiK{\pi) p{pi) ^D{X) , (2.16) gsie 

J-A 



where the integral kernel K and the driving term D are given by 

r^f sm{iii) -sin(2//) 



sinh(// — A -|- 2iri) sinh(// — A — 2irj) ' sinh(A — irj) sinh(A -|- irj) 

(2.17) kernel 

For the XXX case we have 

2 1 
^xxx(/^,A) = - _ ^ 1 ' ^xxx(A) = -j^^ ■ (2.18) kernelXXX 

Here A depends on the external magnetic field h. The physical picture of the ground state 
is that of a filled Fermi sea with boundaries ±A. The dressed energy of a particle in the sea 
is given by the solution of the integral equation 

e(A) t dn K{X, ii) e(/x) =2h- ^^^.^y.^^^^^' — . (2.19) gsie2 

^ ' 2nJ_j^ ^ ^ sinh(A - i?7) sinh(A + i77) ^ ^ ^ 

The requirement of the vanishing of the dressed energy at the Fermi boundary e(±A) = 
determines the dependence of A on h. For small h this relation can be found explicitly by 
means of a Wiener- Hopf analysis '^^1 . For h > — (2 cos 77)^ the system is in the saturated 
ferromagnetic state, which corresponds to A = 0. 
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3. Two-Site Generalized Model 



For the evaluation of correlation functions the so-called "two-site generalized model" 
has proven an extremely useful tool. From the mathematical point of view this is simply the 
application of the co-product associated with the algebra defined by (2.7). The main idea is 
to divide the chain of length L into two parts and associate a monodromy matrix with both 
sub-chains, i.e. 

= T(2,/.)T(1,/.) , T(.,/.) = (^^jII^I ^j^l) {i = 1,2) . (3.1) tl2 

In terms of L-operators the monodromy matrices are given by 



T(2, /i) = L^{fi)L^_^{fi) . . . L„(/i) , 
r(l,//) =L„_i(//)L„_2 (//)... Li(//) . 



(3.2) tl21 



By construction it is clear that both monodromy matrices T{i, /i) fulfill the same intertwining 
relation (2.7) as the complete monodromy matrix T[ii). Similarly the reference state for the 
complete chain is decomposed into a direct product of reference states |0)j for the two sub- 
chains |0) = 10)2 (E) |0)i. The resulting structure can be summarized as 

A(/.)|0), = a,(/x)|0) , A(/^)|0). = rf.(/^)|0), , 

(3.3) abcdi 

C(/.)|0), = 0, S,(/x)|0),^0, 

where the eigenvalues a and d in (2.12) are given by a(/^) = a2{ij)a^{ij) and d{ii) — 
^2(a*)^i(a*)- The creation operators -B(/x) for the complete chain are decomposed as -B(//) = 
^2 (a*) ® -^i(a*) + -^2 (a*) ® -^i(a*)) which implies that eigenstates of the transfer matrix can 
be represented as 

N "i nj 

nB(A,)io) = En n <^i^'i)<iMV) _ 

j=i i,ii jei ken v"^-^^ 

X /(Aj, A^O (^2(Af )|0)2) ® (i?i(Aj)|O)0 , 

where the sum is over all partitions {Aj} U {A^"^} of the set {A^} with card{A^} = n^, 
card{A^^} = n2 — N — n-^. A similar equation holds for dual states 

N ^2 

j=i i,ii jei ken V-^-^^ 

x/(Af,Aj)G(0|C,(Aj))®(2(0|C2(Af)) . 
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4. Reduction of Correlators to Scalar Products 



In this section we reduce the problem of evaluating correlators of the form (cr|(j^) (where 
denotes the normalized zero temperature vacuum expectation value, i.e. the expectation 
value with respect to the antiferromagnetic ground state described by (2.16)-(2.19)) to the 
computation of certain scalar products between states given by the Algebraic Bethe Ansatz. 
We start by noting that due to translational invariance it is sufficient to consider the cor- 
relator G{m) = (c^Ci). In terms of the operators Qj = |(1 — crj) the correlator takes the 
form 

G{m) = ^q^q,) - A{q,) + 1 , (4.1) g 

where we have again used translational invariance. The quantity (g^) is nothing but the 
density of down spins in the ground state and can thus be reexpressed as {q^) = /^^ dAp{A). 
The first term in (4.1) is expressed in terms of the quantity Qi{m) — Yl"Li Qj = Sjli '^J'^t 
as follows 

(QmQi) = lAiiQAm))') , 

where A (not to be confused with the inhomogeneity A in the XXZ hamiltonian) is the 
lattice laplacian A/(j) = /(j) + /(j — 2) — 2/(j — 1). Putting everything together we obtain 

r-A 



G{m) = 2A{{Q^{m) f) + 1 - 4 / dX p(A) . (4.2) szsz 

J-A 



The only nontrivial quantity to determine is thus {{Qi{'m)y) — {exp{aQ-^{m))) 
will now use the two-site generalized model to express the "generating functional" 



. We 

a=0 



„ . , , n, ^^^ <"! H,".! C(Aj) exp(a Er.i -^r-^D ULi , 
:^ (exp(„«.(™))) := (0| ^ . C(A,) i^(AJ|0) ' ^''^ 

in terms of scalar products: we take the first sub-chain to contain sites 1 to m and the second 
one sites m -|- 1 to L. We note that Qi{m) now acts only on the first sub-chain and simply 
counts the number of down spins. Using (3.4) and (3.5) in (4.3) we obtain 



F(a, m) = i(0| n C,{X%) [J B,(XfJ\0), ,(0| J] C,{Xf,J [J ^2(Af,J |0), 

^ Ic Ib lie IIb 

X n «2(Afjrf2(A?j n «i(A?/j^i(Af,j n /(^f^'^f/j n /(aF/,,a?j, 

IbjIc IIbJIc IbJIb Ic'IIc 

(4.4) genfu2 

where the sum is over all partitions 
{AfJ u {Xf,J = {A}, {AfJ n {Af,J = , {Xfj U {A,%} = {A}, {x^} n {Xf,J = 
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of the set {A} with card{Aj^} = card{Aj^} = n^, card{Ajj^} = card{Ajj^} — N — n-^ and 

N N 

<^iv = (0|n^(\)n^(^^)|0) ■ (4-5) Sn 

j=i k=i 

Note that due to (3.3) and (2.7) (for B-{fj.). C^{fi)) wc only need to consider partitions such 
that the size of partitions and Iq< (and 11^ and He) arc the same. In the following 
section we will show that scalar products of the form appearing in (4.4) can be expressed as 
determinants and then use this fact to obtain a determinant representation for F(Q;,m). 

A particularly simple correlator to compute within this approach is the "Emptiness 
Formation Probabihty" . This correlation function is defined as 

and physically corresponds to the probability to find a string of ferromagnetically ordered 
adjacent spins in the (antiferromagnetic) ground state'^^l. It can be obtained from F{a^m) 
in the limit o; — > — oo 

P{m) — lim F{a,m) . 

a-^—oo 



5. Scalar Products 

We now turn to the investigation of scalar products of the form 

N N 



SN-{0\llC{Xf)U^iX!m . (5.1) scpr 



fe=i 



Here we do not assume that the sets of spectral parameters {A"^} and {A*"} are the same, 
and we also do not impose the Bethe equations (2.15), because out goal is to determine the 
scalar products occurring in (4.4). The norms^^^ are a special case of (5.1). Prom (2.8) and 
(2.12) it follows that scalar products can be represented as 



= E n <^f) n di^^)^N ( l^^l l^^l ) , (5.2) Ks 

where the sum is over all partitions of {A'"} U {A^} into two sets {A^} and {A^}. The 
coefficients iCj^ are functions of the A^ and are completely determined by the intertwiner (2.7). 
In particular the -f^^v's are identical for the homogeneous model (2.6) and the inhomogeneous 
model (2.11), i.e. the Kj^''s are independent of the inhomogeneities {u^} and also do not 
depend on the lattice length L as long as N < L. We will exploit this fact by considering 
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special inhomogeneous models for which all terms but one in the sum in (5.2) vanish, and 
then represent this term as a determinant. The basic tool for representing scalar products 
as determinants is a Theorem due to Izergin, Coker and Korepin'^*'^^^, which deals with 
determinant representations for the partition functions of inhomogeneous XXZ and XXX 
models constructed according to (2.11): 

Theorem 1: Consider an inhomogeneous XXZ chain of even length N with inhomogeneities 
j = 1...N. Let |0) and |0) be the ferromagnetic reference states with all spins up 
and down respectively. Let B{fi) and C{fi) be the creation/ annihilation operators over the 
reference state |0) . Then the following determinant representations hold for the XXZ magnet 

N N 

(0|n^(A,)|0) = (0|nC(A,)|0) 

TV TV 

= (-1)"^ n n ^^"^HK -^k- in) sinh(A, + ir)) (5.3) detZl 

a=l k=l 

X Yl sinh(A„-A^) Yl sinh(z/; - z/J j det{M) , 

\l<a<f3<N l<k<l<N / 



where 



. > i sin(2?7) , ^ ^, 

^ak = — TT^ . , . ' ■ 5.4 detZ2 

smh(A^ - z/fc - tr]) smh(A^ -i'k + W 



A similar representation holds for the XXX magnet. 

Let us now derive explicit expressions for the coefficients Kj^. It will be convenient to 
work with the following sets of spectral parameters 

{A^^} = {A^} n {A^}, {A^^} = {A^} n {A^}, 
{A^^} = {A^} n {A^}, {A^^} = {A^} n {A^}, 

with cardinalities 

n = card{A^^} = card{A^-^} , 

N -n = cardlA-^^} = card{A^^} . 

The partition with n = is characterized by {A^^} = {A^}, {A^^} = {A^}, {A^-^} = = 

{A^*-^}. The corresponding coefficient Kj^ ( r\C'i r\Bi ) called highest coefficient. 

V|A I |A 1/ 



Lemma 1: For highest coefficients the following determinant representation holds 

(a^l an) = (n^(^-^")^(^^'^'')) nMA?,Af)det(M-) 



/ ^■''^ (5.5) Kcbcb 

h(u lA - iM''\ - - f(\c xB) 
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For the XXZ magnet we find 

h(X,Li) = ^ f-- ^ , t(X,n) = — — ^ ^ [" , r . 5.6 Kcbcbxxz 

^ '^^ ism{2ri) ' v >py sinh(A - + 2ir/) sinh(A - /x) ^ ^ 

and in the XXX case we have instead 

/i(A,/x) = 1 -i(A-/i) , t(A,/i) = -— --^ — . (5.7) Kcbcbxxx 

Proof: We will carry out the proof for the XXZ case, the XXX case is similar. Consider 
an inhomogeneous XXZ model on a lattice of length N with inhomogeneities = A^ + ir]. 
We have a(A) = H^i sinh(A — Xf — 2irj) and d{X) = H^i sinh(A — A^). Inspection of (5.2) 
yields that in this situation only one term in the sum of the r.h.s of (5.2) survives, namely 
the one with {A^} = {A^}. Thus for this special scalar product we obtain 



S 



N 



Uj=Xf+ir] 



( \P\ I An ) n «i^h(Af - A^ - 2trj) n sinh(A^ - Xf) . (5.8) hcoef f 



On the other hand -B(A) flips one spin, and as we have chosen N to be the length of the 
lattice we find that IljLi -^(-^jOlO) proportional to the ferromagnetic state with all spins 
flipped, and thus orthogonal to all states in a basis other than |0). Thus 



N N 

Sn ={0\UC{Xf)\0){0\l[B{X^m . 

i^j=Xf+iV j=l k=l 



By Theorem 1 both factors can be represented as determinants. By direct computation we 
find for one of the factors 

N 

(0| n^^(Af )|0) = n«inh(A^ - xf - 2vn) . 

Using the determinant representation given by Theorem 1 on the other factor we arrive at 
(5.5). q.e.d. 



Lemma 2: Aihitiaiy coefficients Kj^ are expressed in terms of highest coefficients as follows 

\\/^ J -j^A J-/ \jeACk€DC / \l<=ABmeDB 

"V{A^^} {A^^i; ^""V{A^^} {A^^}y 



(5.9) Kcbad 



Proof: We again will only treat the XXZ case explicitly, the XXX case being very similar. 
Consider an inhomogeneous XXZ model with inhomogeneities {uj} — {Xf^+iriyulXf^+irj}. 
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Now only the term proportional to Kj^ 1 {A"^} {A^}J ^^^'^ 
survives. Proceeding as in the proof of Lemma 1 above we arrive at (5.9). q.e.d. 

Combining the results of Lemmas 1 and 2 with (5.2) we arrive at the following expression 
for general scalar products of XXZ and XXX magnets 

^iv = n ^(^''' ^^)^(^^' ^f) E «gn(^c)sgn(PB) n MAf , ) n MA^ , AD 

i>'^ i'^ ^."1 Sn2 

X nMAr'^^'')nMAf,Ar)det(M^^)det(M^g), 

l,k j,m 

where P^* ^^e permutation {A^^*^, . . . , A;^*^, Af^, . . . , A^'l„} of {Af,...,A^}, is the 
permutation {Af ^, . . . , X^^ , A^^, . . . , A^^„} of {Af , . . . , A^}, sgn{P) is the sign of the per- 
mutation P, and 

(^^c),. = t(Af , Ar)rf(Ar)a(Af ), t(A,/.) = i^;^ . 

Note that (5.10) is formally the same as the corresponding expression for the delta-function 
Bose gas^^^l, only the functions /(A, /x), g{X,fi) (and thus also h and t), a(A) and d{X) are 
different. 



6. Dual Fields 



The most important step in the DFA follows next: we introduce dual quantum Gelds 
in order to simplify (5.10) and obtain a manageable expression for scalar products. This 
step was first carried out in [6] for the delta-function Bose gas. The XXX and XXZ cases 
of interest here can be treated very similarly, so that we will be brief in our discussion. The 
fundamental observation is that the r.h.s. in(5.10) looks like the determinant of the sum of 
two matrices: 

Lemma 3: Let A and B be two N x N matrices over (D. Then the determinant of their 
sum can be decomposed as follows 

det(^ + E) = J^sgn(Pjsgn(Pjdct(y4p^pJdct(Ep^pJ . (6.1) detAB 



Here P^ and P^ are partitions of the N rows and columns into two subsets TZ, TZ and C C 
of cardinalities n (for TZ, C) and N — n (for TZ, C) respectively, Ap p is the n x n matrix 
obtained from A by removing all TZ-rows and C-columns, and Bp p is the N — n x N — n 
matrix obtained from B by removing all TZ-rows and C-columns. Finally sgn(Py,) is the parity 
of the permutation obtained from {1, . . . ,N) by moving all TZ-rows to the front. 



11 



Proof: See [14] p.221 ff. 



Comparison of (5.10) with Lemma 3 shows that one does not get the h{X, /i)-factors by 
simply taking the determinant of the sum of the matrices M^-^. This leads to the introduction 
of two dual quantum fields $yi(A) and $£)(A) which are represented as sums of "momenta" 
P4 and "coordinates" as follows 

^aW = QaW + PdW, $z5(A)=g^(A) + P^(A), (6 2) DQF 

[PoW,Qoif^)]^Hh{\f^)), [P^(A),g^(/.)] = ln(M//,A)) . 

All other commutators of P's and Q's vanish. A very important property of the fields $ is 
that they commute for different values of spectral parameters 

[<I>^(A), = = [$a(A), $a(/^)] = [<f d(A), ^nif^)] ■ 

The dual quantum fields act on a bosonic Fock space with reference states |0) and (0| defined 
by 

P,(A)|0) = 0, (0|Q,(A)=0, a = A^,(0|0) = l. (6.3) dfs 

Using the dual fields it is now possible to recast (5.10) as a determinant of the sum of two 
matrices 

Theorem 2: Scalar products for ttie Heisenberg XXZ and XXX magnets can be represented 
as determinants in tfie following way 

Sn = n ^(^?' Afe )^(Af , Af )(0| det S\0) , 

S,, = t{Xf, Af )a(Af )d(Af ) exp ($^(Af ) + $^(Af )) ^^.4) Sn3 

+ t{X^, Xf)d{Xf)a{X^) exp ($^(Af ) + $^(Af )) . 
Proof: Using Lemma 3 to expand the determinant in (6.4)we arrive at 
(0|det5|0) = ^sgn(P^)sgn(P^)det(M^^)det(M^g) 

(n N-n \ (6.5) pi 

^ $^(Af ) + $^(Af ^) + J2 ^a(A^^) + ^ni^n |0) . 
j=l k=l / 

Evaluating the expectation value of the dual quantum fields by means of (6.2) and (6.3) we 
arrive at (5.10). q.e.d. 

It is possible to further simplify (6.4) by eliminating one dual field: we define a new 
dual vacuum (0| according to 

(0| = (0| exp (j2 PDi>^f ) + ^A(Af ) j , (0|0) = 1 , (6.6) ndv 
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and a new dual field 

<p{x)^p{x) + q{x), q{x) = g^(A) - g^(A) - (o|g^(A) - g^(A)|o), 

p{X) = PM - P^iX) , (0|g(A) = =p(A)|0) , 

[p(A),g(/x)] = -ln(/i(A,/x)M/^,A)) , [p(A),p(/x)] = = [g(A), , [99(A), (^(/x)] = . 

(6.7) phi 

In terms of this field we obtain the following determinant representation 

N 

Sn = n^?(\'''A^)^(A^\'')n«(^'')^(^'')nMAf,Af)(0|det5|0) , 

j>k j=l j,k 

S,k = t{Xf, Af ) + t{X^, Af )^ exp ((^(Af ) - ip{Xf)) (6.8) SnSnew 

^(Af , ^m)H^my A?) 



X 

m 



n 



where r(A) = Now we have all the machinery ready to tackle the problem of representing 
(4.4) as a determinant. 

7. On Norms 

In this section we will have a closer look at norms of Bethe wave functions. These were 
evaluated for both XXZ and XXX models in [5], so that the answers are already known. Here 
we will consider norms as special cases of scalar products in order to build up some machinery 
needed below for further analysis of (6.4). We will treat the XXZ case in detail and quote 
the results for XXX. In order to study norms we ought to set {A*"} = {A^} in (6.4) and 
then impose the Bethe equations (2.15). Immediately some problems arise as the diagonal 
elements of the matrix S in (6.4) become ill-defined ("^") and have to be investigated more 
carefully. The off-diagonal matrix elements are easily dealt with. The Bethe equations (2.15) 
together with the antisymmetry property g{X, fj,) — —g{iJ>, A) imply 

^(^^) n WTyI - (-l)""" ,k^l,...,N. (7.1) bae2 

Thus we obtain 

S^j, = i(A^., A,) + t(A„ A^.) exp (<^(A,) - ^{X^)) , j k . (7.2) od 

To obtain the diagonal matrix elements we take the limit A^- — > Xj^ in the matrix S and use 
I'Hospital's rule (here we have to make use of the explicit expressions for the functions /, g, 
a, d, etc for the XXZ case) 



d 

Sjj = zsm{2r]) — 



^(A, A„^ ^ 



dip 

— 2 cos(2?7) + i sin(2?7) — 

A=A, ^A 



. (7.3) d 
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To obtain this expression we also have made use of the Bethe equations (7.1). We observe 
that the last two terms in (7.3) are precisely what one obtains when taking the limit — > 
in (7.2). Putting everything together we find the following expression for the norm (4.5) 

N N N 

j=l k=l jj^k j=l 

J\fjk = t{Xj, Xk) + t(Afc, A^.) cxp (v9(Afc) - <^(A^.)) (7.4) norm2 



d 



ln(r(A)) + J]ln( 



n=l 



MA„,A)' 



A=A, 



where we now interpret the first two terms in Af in the sense of I'Hospital for the diagonal 
elements. For the case of the XXX magnet we have to replace sin (277) by 1 and use the 
functions f,g,h,t following from (2.3). There is one further simplification: it was shown in 
[5] that the expectation value of the dual field part in (7.4) is such that the dual fields can 
be simply set equal to zero, i.e. we can replace (0|detA/'|0) by dctA/"', where J\f' is obtained 
from jV by dropping the exp((^)-terms. Then a further simplification takes place as 

t{X^,X,)+tiX„X^) = -sm{2r])K{X^,X,), 

where K is defined in (2.17). This is summarized in the following theorem due to Korepin'^l 



Theorem 3: Norms for the Heisenberg XXZ and XXX magnets can be represented as 
determinants in tlic following way 



N N N 

j=l k=l jj^k j=l 



For the XXZ case the matrix M' is given by 



AT' = sin(277) -i^(A,, A,) + ^ 5 



_d_ 
'^''dX, 



N 



n=l \ rt' . 



MA,, A J 
■) 



sm{2rj)[-K{X^,X,) + 6j, 



_d_ 
dX^ 



N 



MKA,)) + $^i^(A,,AJ 



n=l 



(7.5) normS 



(7.6) normSxxz 



where K{X, ji) and h{X, /i) are defined in (2.17) and (5.6) respectively. For the XXX case we 
have instead 



jk 



(A, - A,)2 + 1 



+ 6. 



jk 



(7.7) normSxxx 



Proof: see [5] . 
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8. Correlators on the Finite Chain 



Let us now come back to the generating functional for correlators (4.4). We will now use 
the machinery built up in the last few sections to express -F(a, m) as a determinant. We will 
proceed in two steps: we first will analyse (4.4) without using that {A^} = {A*"} = {A} and 
without imposing the Bethe-equations (2.15). In the second step we will then impose these 
two constraints. Using (6.8) we can represent the scalar products in the two-site generalized 
models in (4.4) as determinants 



Lemma 4: 



n MAF„A?)(0|det.,({Af},{Af})dets,({AF,},{Af,})|0) , 

n N—n 



{s,{{X% {A^})),, = exp(a)d,(Af )a,(Af ) (.".({A^}, {A^})) , 
is,i{X% {X^}))j, = a,{Xf)dM {s,{{X% {A^}))^., , 

where 



(8.1) genfuS 



(8.2) genfu4 



(^({A^}, {A^}))^., = t{Xf, X^)a^{Xf)d^{X^) exp ($^^(Af ) + $^^(Af 

+ t{Xl Xf)d^{Xf)a^{X^) exp ($^^(Af ) + $^^(Af )) 

Here the dual fields are defined according to 

^A,W = Qa,W + Pd,W, $D,(A) = Qd,(A) + Pa,(A), 
[Pd,W,Qd,{i^)] = S^^\n{h{X,i,)), [Pa^{X),Q^^{i,)] = 5,^ ln(/i(/., A)) . 

All other commutators vanish. The reference state |0) and its dual (0| are annihilated by all 
momenta/coordinates respectively 

P„(A)|0) = 0, (0|gjA) = 0, a = A^,Z^^,(0|0) = l. (8.4) df2 



(8.3) df 



Proof: We use (6.4) to express both ^{0\ H/^ C^{XfJ Ui^ B^i>^fj\0)^ 

and 2(01 ri//^ ^2(AF/c) ri/Za -^2(Af/^)|0)2 as determinants. Wc arc led to introduce two 
sets of dual fields (one for each scalar product) $^ (A), $^ (A), 7 = 1,2 with commutation 
relations given by (8.3). The two kinds of dual fields are completely independent of each other 
(all commutators bewteen momenta/coordinates of different sets vanish). The representation 
(8.1)-(8.2) is now obtained by direct computation, where the sgn(Pg)sgn(Pf;) arises upon 
taking the factor nj>fe fl'(Aj^, A^)gr(Af , A^) in front of the sum due to g{X,iJi) — —g{ijL,X). 
q.e.d. 
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We now observe that (8.1) is basically of the same structure as (5.10). Thus, in analogy 
with (6.4), we can introduce new dual quantum fields and reexpress F{a,m) as a single 
determinant. 



Lemma 5: Consider the set of four commuting dual quantum fields 

^D, (A) = Qd, (A) + P^, (A) , (A) = Q^^ (A) + Vn, (A) , 
^D,W = Qi.,(A) +^A,(A) ,V'A,(A) = QA,W+rn,W , 
with commutation relations of the momenta/coordinates given by 

[^D^W, Qd,{i^)] = 5,^ln(/i(A,/.)), [P^^(A), Q^^(/.)] = 5^^1n(/i(/., A)) . 

All other commutators vanish. The action of the dual fields on the dual reference states is 
given b7 7'„(A)|0) = 0, (0|Q„(A) = 0, a = A^^A^^D^^D^. Then the following determinant 
representation holds 

F{a, m) = — n 5(Af , A^)^(Af , Af ) (0| det M |0) , 

+ (^2({A^}, {A^}))^.,exp (^^^(Af ) + ^^^(Af )) , 
where (s^)jfe are given by (8.1). 

Proof: The proof is analogous to the one for Theorem 2, only the expectation value of dual 
quantum fields is slightly different, q.e.d. 

So far wc have not used the fact that we are dealing with expectation values of Bethe 
states, i.e. we have neither used the fact that {A*"} = {A^} = {A} nor imposed the Bethe 
equations (2.15). In the next step we will impose these constraints. The discussion will be 
reminiscent of section 7 above. The result is summarized in the following 

Theorem 4: The generating funtional F{a, m) can be represented as a ratio of determinants 
in the following way 

^ ' > det AT' ' 

r (A ■) 

^ok = ^(Aj, K) + ^(Afe, A^.)^— ^ exp ((^2(Afe) - <^2(Aj)) 

\ ^i(A) 
+ exp {a + (^4(AJ - <^3(Aj.)) \) + ^(Aj, Aj^)^— ^ exp ((^i(Aj.) - (^i(AJ) 



-r5^,u ||LD(A,.)+^X(A,,Ajj, 



(8.6) genfu6 
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where r^{X) = a^{X)/di{X), K{\n) and D{X) are defined in (2.17) (2.18), u = sin(277) for 
XXZ and u; — —1 for XXX, and the commuting dual fields are defined according to 

VaW = PaW + QaW , i^QaW = = p„(A)|0) , (0 10) = 1 , a = 1 . . . 4 , 

/I 1 0\ /I iN^ 



10 1 
111 

Vi 1 i/ 



ln(MA,/x)) + 



110 
10 11 

Vo 1 1 i/ 



ln(/i(/x, A)), a,6 = 1 . . .4 . 



(8.7) ndfS 

Here all terms not proportional to Sj,^ in are understood in the sense of I'Hospital for the 
diagonal elements. 

Proof: We start by defining a new dual vacuum (0| and a new set of dual fields according 
to 



N 



, (0|0) 



(0| = (0| exp PoS^i) + PaS^^) + 'PoSh) + ^A,(A 
\j=i 

<t>,{X) = $^^(A) - $^^(A) , UX) = $^^(A) - $^^(A) , 

U>^) = ^A,(A) - ^D,(A) - + $^^(A) , 

04(A) = V'A,(A) - ^z.,(A) + <I>^^(A) - $^^(A) . 

The fields (paW can be decomposed into momenta and coordinates g„ (by using (8.3) 
and the definitions of V'a given in Lemma 5), which are found to obey the commutation 
relations (8.7). By straighforward rewriting of (8.5) in terms of the new fields and the new 
dual reference state we obtain 



8.8) ndv2 



.9) ndf 



N 



F{a,m) = — n/(A,-,A,)/(A„ A,.)n«(A,)^^(A,)(0|detA^|0) , 
'^^ j>k 

M^,^tiXj,X,)+t{X„xf^ 



exp (02(Afe) - (/)2(Aj)) 



+ 



KA,) 



exp [a + 04(Afe) - 03(A^)) 



^(Afe, A ) + t{Xj, Aj._ 



^i(A.) 
^i(Afc) 



exp(0i(A^.) -0i(Afe)) 



(8.10) genfu7 



Here we have used that 



N 



(0| exp <t>A, (A,) + <t>D, (A,) + (A,) + ^n, (A,-) = \{ h{Xj, AJ (0| 



It is found that whereas Pa{X)\0) = 0, the coordinates ^^(A) of 0(j(A) do not annihilate 
the new dual reference state {0\. Therefore we "shift" (p^W by subtracting their vacuum 
expectation values in analogy with (6.7) 



faW = <PaW - (O|0„(A)|O) =p„(A) +g„(A) , a = 1, 
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(8.11) ndf2 



By construction the p's and q''s have the same commutations relations (8.7) as the mo- 
menta/coordinates p^(A) and of the 0„(A)'s. Furthermore p„ ( A) |0) = and (0|g„(A) = 
for a = 1 ... 4. The shifts are found to be 



«:„(A) = (O|0,(A)|O) = (l-<5„,)5]ln 



^(A, A^.) 
MA., A) 



If we replace the fields 0^ in (8.10) by the fields (p^ we pick up additional factors due to the 
shifts 

1 ^ 
F{a,m) ^ —llf{X^,X,)lla{Xj)d{X^){b\detg\0) , 

g.^ = t(A., AJ + i(A„ A.)^^e'^2(^'=)-'^2(^^)e'^2(^'=)-'^2(^^) 



r(A,) 



^(Ajfc, A) +t(A-, Afc 



ri(A, 



exp (<^i(A^) -(/?i(Afe)) 



^ ^' '^^^(A,) 

(8.12) genfuS 

The off-diagonal matrix elements of G can be further simplified by simply imposing the 
Bethe equations. Rewriting the Bethe equations (2.15) as 



.(A,) n 



3 

N 



MAfe,A,) _ (-1)^-1 
\h{X^,X^?) ri(Afe) 



^.13) baeS 



'^2(Afe) h{Xf^,Xj) 



-1)--V,(A,),^ = 1, 



TV 



we find that the additional factors take the form 



exp(«;2(Afe) -«;2(Aj)) 
exp {nMk) - i^ii^j)) 



r^{X^) ri(Afe) 
KAJ _ . 
KA,) ■ 



Inserting this into (8.10) we arrive at (8.6) without the term proportional to 5^^., i.e. we have 
proved (8.6) for the off-diagonal matrix elements. To get the diagonal matrix elements we 
have to investigate the limit A^- — > A^. of (8.12) in detail. In the limit A^ Xj. the sum of 
the first two terms in and the expression in brackets are both of the form By using 
FHospital's rule we find analogously to section 7 above 



lim U(A^.,A,)+t(A„A^.) 
= —2 cosh(2i?7) + sinh(2i?7) 



^2(Afc) 

'''r,{X,) 
5v^2(A) 



dX 



sinh(2i?7) 



_d_ 
dX, 



Mr2(A,.)) + EM 



^(A,,AJ 

HK, A,) 
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= -2cosh(2,,) + =i„h(2,,) ( ^^H^ + ""'wt, ■ 

^ ^ \ dXj smh{Xj + IT)) sinh(Aj - irj) J 

Using these expressions we find that the diagonal terms of G in (8.12) are equal to the 
diagonal terms of G in (8.6) if we keep in mind that the first two lines of Gjf^ in (8.6) are 
interpreted a la I'Hospital for j = k. Last but not least we insert the expression (7.5) for 
0"^ in the resulting expression and arrive at (8.6). This completes the proof of the theorem, 
q.e.d. 

Theorem 4 states the determinant representation for the generating functional F{a, m) 
on a Unite chain of length L. As always in Bethe Ansatz solvable models significant simpli- 
fications take place if we take the thermodynamic limit L — > oo. This is done in the next 
section. 



9. Thermodynamic Limit 



The results of taking the thermodynamic limit of (8.6) and main results of this paper 
are summarized in the following 

Theorem 5: In the thermodynamic limit the generating funtional F{a,m) for the case of 
the XXZ magnet can be represented as a ratio of determinants of Fredholm integral operators 
(^id + and (^id — in the following way 

{0\det{td +^V^)|0) 
F{a,m) = — — . (9.1) genfu9 

det[zd — 2^-f^) 

Here (0| and |0) are the vacua of the dual bosonic Fock space defined in (8.7) and the integral 
operators act on functions f defined on the interval [—A, A] according to 



{id - ^k) * / 
(id + i-\/)*/ 



where the kernel K{X,ii) is defined in (2.17) and the kernel ofV is given by 
V{Xii)- ~^^^^2r?) / 1 e2^(A)e2(/x) 



sinh(A — /i) 1^ sinh(A — n + 2ir)) sinh(/x — A + 2i7]) 
+ exp(a + ^,(/.) - ^3(A)) ( . -\ , ^. , + ^"(^)^^(^) 



sinh(// — A + 2ir]) sinh(A — + 2ir]) 



(9.2) V 
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The dual fields (p^W are defined in (8.7), with /i(A,/i) given in (5.6). 

Proof: We begin by taking the thermodynamic hmit for the norm cTjy (7.5). We first write 
J\f' as the product of two matrices: 

K 

jk — Sin(2?7) Ijm'^mkJ ^ jm ~ ^jm 'a ' ^ jm ~ ^jrrfirm 

m 

where 6^ = LD{\J + E„^(^m>\)- Here D and K are defined in (2.16)-(2.17). The 
determinant of A/"' is the product of the determinants of / and J. Next we use that the set 
of roots {Xj} describes the ground state and the roots thus obey the equations 

N 

27rLp{X.) - ^ K{X^, A,) = LD{X^) , j = l...N , 

k=l 

which is the discrete version of (2.16). Here p(A •) = Yn—^—ri- which becomes p(A) defined 
by (2.16) in the thermodynamic hmit. We thus can rewrite 6^^ = 27rLp(A^), which leads to 

TV 

det J = 27rLp{Xj) . (9.4) J 

i=i 

In the thermodynamic hmit the matrix / turns into an integral operator I — id — 



2ir 

A 



A 



f{X)-^J^^di,K{X,p)f{i,), 



where K is the kernel of K defined by (2.17). 

The matrix G^-^ in (8.6) is treated in a very similar way. We rewrite it as a product 

Gjk = sin(27;) J2 ^jmJmk > 

m 

where = 6j^27TLp{X^) is the same as above, and 

r ^i('^7) \i 1 

+ exp [a + (^4(Afc) - <^3(A^.)) ^{X^, Xj) + t{Xj, A^)^;-^ exp {^i{Xj) - ViiK))^ j ■ 

In the thermodynamic limit the matrix Wji. turns into an integral operator W — id+ -^V, 
with kernel y(A, //) defined by (9.2). Thus we obtain (9.1) in the thermodynamic limit, 
q.e.d. 
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For the XXX chain a determinant representation is obtained in an analogous way. The 
result is found to be 

Fxxx{»,m) = -— ^ — , (9.5) genfuxxx 

det(l - ^Kxxx) 

where -f^xxx ^xxx integral operators with kernels K{X,ii) from (2.18) and 

A — /ilA — u + z Li — A + t , , 

/ -1 / N / X X ^ (9-6) Vxxx 

.exp(...,.)-.3(A))(^.^)} 

exp((^,(A)), e,(A)=(^^j exp((^,(A)) . 

The dual fields ^^^(A) are again defined in (8.7), but now h{X, n) = 1 — i{X — jj). 

The Emptiness Formation Probability can be easily obtained from (9.1) and (9.5) by 
setting a — —oo, which corresponds to dropping the second line in the expressions for the 
kernel of V in (9.2) and (9.6). For the XXX case this exactly reproduces the result of [4]. 



10. Some limiting Cases 



It is quite straightforward to evaluate the determinants in (9.1) for strong magnetic 
fields h ^ hf, — (2 cos 77)^, in which the ground state is very close to the ferromagnetic 
vacuum and A <S 1. The near asymptotics (m <^ (7r/2A) tanry) of the EFP for the XXZ 
case follows to be 

P(m) = 1 — I cos 77 I m . 

\7rsinr7 / 

Using (2.16) and (2.19) this reproduces the obvious result 



P(m) 



2 ^ 



1 — —\/h^ — h , h 

TT 



Another interesting limiting case (which allows to make contact with known results) is 
the XXQ free fermionic limit of the XXZ model, where = ^ in (9.2) and (2.17). This case 
has been previously considered in [20,21], where a determinant representation for F[a,m) 
was found, which does not involve dual quantum fields. Taking the free fermionic limit of 
(8.7) we obtain 



[%(/^),Pa('^)] 



/2 


1 

Vi 




2 
1 
1 



1 
1 
2 
2 



1 
2 

2/ 



In (cosh (A — /jl)) . 
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Thus we can choose (p^{X) — (f^iX) and reduce the number of dual fields to 3. Furthermore 
we have sin(477) = and therefore K{X, ji) — 0. The determinant formula then reads 

F;,^(a,m) = (0|det(io?+^T/)|0), 

^^^^-^""^ ^ sinh2(A-/.) r " '^"'^^^'^^^^ " e«+^3{.)-^3(A) ^ _ e-,\lAe,{X)^ | , 



e —i\ , ,,,, (e +1^ 



where now 

^2(^) = (^^^A^ J exp((^2(A)) , ei(A) = (^^^x^i J ^MVi{X)) ■ 

This expression has to be compared with the result obtained in [20] which reads after trans- 
forming their expression to the notation used in the present paper 

F-^-^{a,m) — det [id -\ U] . 

Here U is the integral operator defined in terms of the kernel 

u(x,,) . ^ilLzf!). U _ (^^, . (10.2) „ 

^ ' smh(A - ^Ji)\ \e^^-i e^i' + ij j ^ ' 

The proof of the equivalence of (10.2) and (10.1) is quite tedious and might be presented 
elsewhere. 



11. Discussion 

In this paper we have derived determinant representations for spin-spin correlation 
functions in XXZ and XXX Heisenberg models. Other correlation functions can be treated 
similarly. In a following publication we will use the determinant representation to embed 
the correlation functions in systems of integrable integro-difference equations^^^l . 

Finally we would like to emphasize that the Dual Field Method can be applied to any 
correlation function in any integrable model. 
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